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ABSTRACT 

X-rays from Cygnus X-3 have been observed during early 1978 with 
Che detectors of the SAS-3 satellite. These observations, in conjunc- 
tion with earlier URURU and ANS data. Indicate that the 4.8 hr period 
of Cygnus X-3 is increasing at a rate of P/P - (5.1 1 1.3) x 10 ** yr * . 
The sign and magnitude for this change are incompatible with a rotation 
model for the period and are in reasonable agreement with model pre- 
dictions for orbital changes associated with mass loss and transfer in 
a binary system. 
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Search for X-Ray Emission from the y-kay 
Source y!95+5 . R. C. IAMB, Iowa State Univ. and 
R. C. DOWER, M. I.T. — The high energy (> 100 MeV) 
y-ray source, yl95+5, observed by both the SAS-2 and 
COS-B satellites with an intensity approximately equal 
to that of the Crab, has so far elluded detection at 
other wavelengths. In an effort to establish a 
possible X-ray identification an observation with the 
SAS-3 rotation modulation collimators was conducted 
January 5-12, 1978. Preliminary results show no X-ray 
source of greater than 6o significance within either 
the SAS-2 or COS-B 90% error boxes, corresponding to 
an upper limit of the X-ray intensity in the 2-11 keV 
range of less than lyJy (''-l Uhuru ct/sec). Comparison 
of this result with unpublished 0S0-8 evidence for an 
X-ray source in this region and recent SAO work will 
also be presented. This work was supported in part 
by the National Aeronautics and Space Administration. 
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I. INTRODUCTIJN 


The 4.8 hour periodic variation in the X-ray int .'natty of Cygnus X-3 

has been observed since late 1970. We present objervatlons from the 

SAS-3 satellite taken during early 1978 which, in conjunction with 

earlier data (Parsignault , et al., 1976, Leach, e£ a \_ . , 1975), show 

evidence for a non-constant period. A simple interpretation of these 

data is that the period of Cygnus X-3 is increasing with time at a rate 
* -6 -1 

P/P ■ (5.1 t 1.3) x 10 yr consistent with the previous authors' 
upper limits. The sign and magnitude for this change are incompatible 
with a rotation model for the period, and are in reasonable agreement 
with predictions by Davidsen and Ostriker (1974) and Pringle (1974) for 
orbital changes associated with mass loss and transfer in a binary 
system. 


II. OBSERVATIONS AND RESULTS 

Cygnus X-3 was observed with the slat detectors of the SAS-3 

satellite (Buff et al., 1977). From 1978 January, 5.9 to 12.0 (UT) the 

source was scanned once per satellite orbit (^94 min.) while another 

region of primary interest was observed with the rotation modulation 

counter system on the satellite's spin axis. During these observations, 

Cygnus X-3 varied from 21° to 27° from the scan plane, and the relative 

efficiency of the center detector varied between 0.39 and 0.32. The 

results of these observations are shown in Figure la, in which the 

source counting rate, corrected for elevation. Is plotted versus the 4.8 

hour phase. Phase 0.0 corresponds to the minimum of the fitting 

function: I * A + A 1 sin [u» ( t - t .)-•?■] as used by Parsignault 

o i rain z 
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et al. (1976), so 1 Is a minimum at t • t n i n * Th* source Intensity ’’as 
determined directly from counting rate plots in which the intensity is 
proportional to the height of a triangle response above a smooth back- 
ground. The peak height varied from ^ 1 to 4 ram above background on the 
plots corresponding to a counting rate of ^ 6 to 24 cts/s. The estimates 
of error Indicated in the figure Include a constant error of + 1/3 mm 
(2 cts/s) in the background determination and a "statistical" error 
proportional to the square root of the peak height. The work of 
Parsignault al^. (1977) has shown that there are intensity fluctua- 
tions larger than those to be expected on the basis of statistics alone. 

Our data qualitatively support this conclusion. In order to achieve a 
2 

X per degree of freedom of *v 1 for the data points near minimum phase 

(O .85 to 0.25) where the excessive fluctuations are minimized, it is 

necessary to increase the "statistical" error by about a factor of two 

from a value which would be consistent with counting statistics alone. 

If this adjustment were not done then the errors associated with fitted 

parameters would be somewhat underestimated. The error estimates shown 

in the figure include this adjustment. 

The values of P = 2ir/u> and t . from the fit restricted to the 

rain 

30 points with phases between 'v- 0.85 and 0.25 are shown in Table 1. 

Values of P and t . from a fit to all 76 points do not differ appreciably 
min 

from these values. The fitted ccjnting rate at phase 0.0 was 20 ± 2 cts/s, 

corresponding to an Intensity over the interval 1.5 to 6 keV of 
-9 -2 -1 

^2.2x10 erg cm s for a Crab-like spectrum. 

The second set of observations span a week beginning 1978 January 
26.8. Four observations of Cygnus X-3 were made near the beginning of 
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Che week and seven near the end. For comparable phases there was no 
significant difference in the intensity. Therefore we have combined 
these 11 observations and plotted them in Figure lb. The errors asso- 
ciated with these points are smaller than the January 6-12 data princi- 
pally because Cygnus X-3 was near the scan plane where detector effi- 
ciencies are nearly 1.0 and the counting rates were correspondingly 
larger. The fitted value of the counting rate at phase 0.0 is 18 i 2 cts/s 
consistent with the early January observations. The values of P and 
C min * rom a *** to ^ points is shown in Table 1. If the fit is 
restricted to those points with phases approximately between 0.85 and 
0.25, then only 5 points remain to determine a function wit^ 4 free 
parameters. This procedure was not considered satisfactory. If A^ 
and are taken from the 11 point fit then a satisfactory fit to P 
and are obtained with only the 5 points between phase 0.85 and 0.25, 

and the values of P and t m ^ n do not differ appreciably from the table 

values. However, the error on t . is somewhat greater than that 

mm 

obtained in the 11 point case, and this error is the one we have chosen 

to list. Both values of t . shown in Table 1 are referenced to the 

min 

barycenter of the solar system; the values of the correction necessary 

are 0*?0024 and 0^0030 respectively. It is assumed that previous 

deterninations of t , include this correction. 

min 

If we assume that the difference in the two values is a whole 
number of periods (106) then a period of 0^19969 + 0.00009 over this 
21 day interval is obtained. This accuracy of 'v 1 part per 2000 is 
comparable to earlier determinations of the period over a similar time 


interval, furthermore the value of the period agrees with previous 
determinations. The near constancy of the period allows one to extrap- 
olate over those intervals of time in which Cygnus X-3 has not been 
observed to determine the total number of cycles from a given epoch 
with confidence that no cycles have been missed nor spurious cycles 
Introduced. Therefore, one can tie together previous observations to 
arrive at a more precise period and test for a period derivative 
(Parsignault et al. , 1976, Leach e^ al . , 1975). We assume the period 
varies linearly with time as: P(t) • P q + Pt. The phase, $, in cycles, 
at an arbitrary time, t, after a reference time, t Q , may be expanded as: 

■ At/P - 4 P(At)^/P^ where At ■ t - t . Our values of t , $, and 
o 2 o o o’ 

the error on $, along with previous determinations, are listed in 
Table 2. Also listed are the values of taken from a least squares 

fit to these observations. The values of the parameters determined 
from the fit are: 

P - 0 < ?1996795 ± 0.0000007 
o 

P/P - (5.1 ± 1.3) x 10" 6 yr' 1 
o 

t o - JD 2,440,949.9173 ± 0.0010 

2 

The residuals from the fit are shown in Figure 2. The x for the fit 

2 _ 

is 12.5 for 14 degrees of freedom. The x for a linear fit with P = 0 

is 29.3. The likelihood ratio of the quadratic to the linear fit is 
2 2 

given by exp[- (X oua j - Xn n )/2] * 4400 to 1, equivalent to a 3.7o 

2 

result on a non-7.ero period derivative. The x for the quadratic fit 
is at the 56% confidence level with none of the data points contributing 
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excessively. This Is evidence chat Che simple variation of Che period 
with time Chat we have assumed Is entirely adequate to fit the seven 
years of observations. In particular there is no evidence for signifi- 
cant discontinuities or jumps in the phase as might be expected if the 
4.8 hr period is the rotational period of a neutron star undergoing 
glitches. 


III. DISCUSSION AND CONCLUSIONS 

The observation of an increase in the period of Cygnus X-3 with a 
time scale P/P : 2 x lO’’ yr places a constraint on various models which 
may be proposed. One possibility which would appear to be ruled out is 
that the 4.8 hour period is the rotational period of a neutron star. 

If the X-rays are from accreting matter from a binary companion onto 
the neutron star then, using the analysis of Rappaport and Joss (1977), 
we estimate that accretion torques would generally accelerate the rotation 
on a time scale of 'V 1 yr. Our observations are clearly Inconsistent 
in magnitude and sign with such a model. 

The most popular view of Cygnus X-3 is that the 4.8 hour period 
is the orbital period of a close binary. If one assumes a circular 
orbit then the separation of the two components is 1.0 x 10^ ^total^Q^^ ^ 
Mass transfer and mass loss will alter the separation and hence the period, 
however there is no unique, model-independent relationship between mass 
changes and period changes (e.g. Thomas 1977; Paczyriski 1971). Of the 
models for Cygnus X-3 in the literature (Davidsen and Ostriker, 1974; 
Pringle, 1974; Basko, Sunyaev and Titarchuk, 1974; Milgrom, 1976; 
van den Heuvel and De Loore, 1973) only the first two give estimates 
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of the time scale for possible period changes. 


In the model of Davldsen and Ostriker (1974) the X-rays originate 
in the accretion of a stellar wind onto a white dwarf, and the wind 
itself modulates the X-rays to give rise to the observed periodicity. 

Their lower limit on P/P is +3 x 10 ** yr * which is in good accord with 
our observation. However, the observation of “v 100 MeV y-rayi from 
Cygnus X-3 (Lamb et^ al. . 1977) has made it more likely that the compact 

object is a neutron star rather than a white dwarf, in view of the 

•• 

evidence that pulsars (neutron stars) can emit y-rays (Ogelman, et al . , 
1976) and the absence of any positive evidence for such emission from 
white dwarfs. If the compact object is a neutron star then its rota- 
tional period may provide an additional modulation to that fraction of 
the X-rays which come directly from the neutron star and are not scattered 
appreciably. Serlemitsos et_ al . (1975) have placed an upper limit of 3Z 
on the fraction of pulsed X-rays with periods in the range of 2 ms to 
several seconds. For periods from 2 to 64 s the upper limit is 10Z 
(Parsignault et_ aK , 1977) however we are unaware of any efforts to 
search for longer periods and suggest that such efforts may be worthwhile. 

The model of Pringle (1974) is similar to that of Davldsen and 
Ostriker in that the primary X-rays result from accretion of a stellar 
wind onto a compact object, and a time scale of 10^ yr Is suggested. 

Thus both stellar wind models accord with the observed value of P/P. 

One feature of stellar wind models is that the observed X-ray 
spectra should be somewhat phase dependent, in particular the absorption 
should vary with the 4.8 hr phase (Davldsen and Ostriker, 1974). 

Previous studies had shown no such effect, however recently Becker et al . 
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(1978) have observed phase dependent spectral changes In the X-ray 
emission which are qualitatively conaistent with such models. 

If we accept the view that the compact object is a neutron star, 

then we find that the X-ray luminosity of Cygnus X-3, at an assumed 

distance of 10 kpc (Lauque*, Lequeux, and Rleu 1972), may be maintained 

-9 -1 

with a mass accretion rate of typically ^3 x 10 M 0 yr . This rate 
is too small by a few orders of magnitude to account for the observed 
period changes. For typical binary parameters the absolute value of the 
fractional change in the period due to mass transfer (accretion) is 
comparable to the fractional mass transferred (Thomas 1977). Therefore, 
if the mass accretion rate we have given is a valid estimate, then one 
can say that the observed changes in the period of Cygnus X-3 are not 
primarily associated with accretion but rather are driven almost exclu- 
sively by mass loss from the system. An estimate of the mass loss rate 
may be obtained by using the relation derived by Davidsen and Ostriker 
for the case where there is negligible mass transfer and where the mass 
lost carries away its angular momentum. Under these circumstances 
M t /M t ■ - y P/P where is the total mass of the system. If we assume 
a system mass of 2 M 0 , then the observed value of P/P implies a mass 
loss rate of ^ x 10 ^ M g yr 

In summary, the observed period increase appears to support a 
binary origin to the 4.8 hr period with reasonable agreement in sign and 
magnitude to the detailed predictions of Davidsen and Ostriker (1974). 

Ve are happy to acknowledge helpful discussions with Saul Rappaport. 
This woik was supported in part by the National Aeronautics and Space 
Administration under contracts NAS5-11450 and NSC 5199. 
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After this work was completed, we became aware of COS-B observations 
(Manzo, Moltenl, and Robba 1978) which, in conjunction with the earlier 
UHURU and ANS data, give a value of fr/P ■ (5.2 + 1.2) x 10 ^ yr . Their 
other parameters, P fi and t Q , are also in excellent agreement with ours. 
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TABLE 1 

SAS-3 OBSERVATIONS OF CYCNUS X-3 
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Time of Observation 

Fitted Period 

Fitted Time of 
Minimum Intensity 
(JD 2.440.00(H) 

1978 January 6-12 

0^19987 t 0.00029 

3517.6420 t 0.0033 

1978 January 26 

0^19947 1 0.00028 

3538.8090 i 0.0090 

-February 2 




SUMMARY OF OBSERVATIONS USED TO DETERMINE P 
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FIGURE CAPTIONS 
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Fi gure la,b : 

The phase dependence of the X-ray intensity from Cygnus X-3 
Jan. 6-12 1978 and Jan. 26-Feb. 2 1978 respectively. Phase 0.0 
corresponds to the minimum of the fitting function: 

1 ■ *0 + A 1 8ln >“ (t - *.!„> - fl- 

Figure 2 : 

Residuals from both a linear fit and a quadratic fit to the obser- 
vations of Cygnus X-3 phase minima by UHURU, ANS, and SAS-3 (thi*» 

work) over a 7 year period. The residuals to the linear fit are 

2 

measured from the dashed line; Xu ■ 29.3. The residuals 

A 1 inear 

from the quadratic fit are measured from the solid curve; 

2 

X . • 12.5. The linear fit assumes P = 0 whereas the quadratic 

fit determines P/P to be (5.1 ± 1.3) x 10 ^ yr 
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